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We present an extensive study on the direct determination of the running coupling a s from the static quark 
antiquark force at short distances, in quenched QCD. We find from our high statistics potential analysis that a q q 
exhibits two-loop asymptotic behaviour for momenta as low as .5 GeV. As a result, we determine the zero flavour 

A rough estimate of full QCD effects leads to the five 



A-parameter to be A^L = .630(38)^ = 293(18)t^ MeV. 



MS 

flavour value OLjjgijnz) = 



102^ A comparison with other lattice results is made. 



1. INTRODUCTION 

Considerable progress has been achieved over 
the past three years in the determination of the 
running coupling from SU(3) gauge theory, i.e. 
quenched QCD on the lattice. Just like in real 
experiment, this is not an easy task, since we 
are dealing with logarithmic effects. Moreover, 
asymptotic scaling is substantially violated in the 
/3-window presently accessible to lattice simula- 
tions. This is to say that the perturbative two- 
loop-formula does not yet relate the lattice spac- 
ing a to the bare coupling, — 6/(3, in terms of 
a constant lattice scale A^: 



aA L 
with 

bo = 



exp 



11 
16^ 



1 



102 
(16tt 2 ) 2 



(1) 



(2) 



Therefore this perturbative relation in the bare 
lattice coupling scheme is not well suited for a 
determination of a s . This scheme is also too far 
off the continuum A/5-scheme to rely on the per- 
turbative recoupling pi 



a L + 5.88a| + C(al), a L = g 2 /4n (3) 



at accessible lattice spacings. 

Three major directions have been followed in 
the past to tackle the problem: (a) abandon the 
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bare lattice coupling for the sake of "improved" 
coupling schemes on the lattice, which are bet- 
ter suited to renormalized perturbative expan- 
sions of short distance operators ||; (b) deter- 
mine ctqq directly from the static gg-force, on suf- 
ficiently large and fine lattices ||f| ; (c) contrive a 
scheme with a volume dependent coupling g(L). 
This novel access to the problem, put forward in 
Ref. ||, is intriguing as it lends itself to the ap- 
plication of finite-size matching techniques. 
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Figure 1. Comparison between the oscillation of 
the potential data (at (3 — 6.4) around V CO nt and 
the expectation of Eq. (8). 



Method (c) has been covered during last year's 
conference with application to the SU(2) case Q. 
The outcome for SU(3) gauge theory is presented 
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Table 1 

The number of independent configurations pro- 
duced and analysed on 16 4 and 32 4 lattices. 



in WOLFF's talk @, while the status of method 
(a) has been reviewed in the contribution of EL- 
KHADRA to this session §]. 

We will focus here on the direct determination 
of a qq from the short distance force and report on 
the status of the analysis of high precision poten- 
tial data in SU(3) gauge theory. In view of the 
quality of today's data, the analysis method |]||] 
is devised to extract ct qq locally at short distances 
rather than deriving it from a global ansatz to 
the force - with short and long distance bound- 
ary conditions - as proposed in Ref. Jl0[ . 

We perform the "classical" Creutz experi- 
ment 0] and extract the potential from the local 
masses 



V(R) = In 



W(R,T) 
W{R,T+l) 



T large 



(4) 



using improved techniques for suppression of ex- 
cited state contributions in order to attain an 
early onset of the asymptotic plateau in T, as 
explained in Refs. [[j^pf]. On more than 2000 in- 
dependent configurations of 32 4 and 16 4 lattices 
(Tab. [I]), measurements of on- and off-axis po- 
tentials have been carried out for 72 and 36 dif- 
ferent qq separations R, respectively, with reso- 
lutions ranging from .03 to .25 fm. UKQCD has 
another set of 50 36 4 configurations at j3 = 6.5 
with mostly on-axis potentials |l5} ] measured. 

The interquark force F is expected to behave 
perturbatively at sufficiently short distances: 



F(R) 



--a qq {R)/R 2 



(5) 



Perturbation theory relates the qq- and the MS- 
schemes by means of a very small first-order co- 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




- B 


° 5 porom. 


z * 


D 4 param. (f=0) 




A l=0 




' 1=1 






— 




: 


t 






^ f 

- 



10 



12 



Figure 2. The sensitivity of the on-axis force at 
/3 = 6.4 on the parametrization of the potential. 

efficient [Ol: 



0«) . (6) 

This connection lets the qq-scheme appear as a 
very promising candidate for a "perfect" lattice 
renormalization scheme 

2. UNFOLDING LATTICE ARTIFACTS 

At present, an a s -analysis of the potential must 
rely on lattice data in the i?-region from 2 to 7. 
This necessitates a careful treatment of lattice ar- 
tifacts. 

One expects [fl5f the continuum and lattice po- 
tentials to be related by a correction term 

V(R) = V cont (R) - SV(R) , (7) 

which should be proportional to the difference be- 
tween the lattice and continuum propagators: 



SV(R) = 7 (AnG L (R) - l/R 



(8) 

As one can see from Fig. |l| (circles), this correc- 
tion exhibits significant oscillations around zero 
for R < 4, with maximal deflection occurring for 
the on-axis lattice potentials. For the determina- 
tion of the parameter 7 in Eq. (||) , we make use of 
two paramctrizations of the continuum potential. 
We follow MICHAEL § and choose as the first 
form 



V cont {R) = V + KR- e/R + f/R 2 



(9) 
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while the second is the Cornell potential from 
Eq. pi (/ = 0). We demonstrate in Fig. [I], that 
this simple ansatz works surprisingly well: with 
the best fit value for j/e = I ~ .6, which inci- 
dentally comes out independently of (3 and /, the 
lattice potential data are seen to bounce around 
the interpolating curve in very much the same 
way as described by Eq. (|§|). 

We emphasize again, that at this stage we are 
merely concerned with unfolding the lattice arti- 
facts, and we better make sure that the unfold- 
ing procedure is sufficiently robust in the sense, 
that the resulting "continuum" force, F, is not 
an artifact of the underlying interpolative ansatz 
Eq. ®. 

To that end, let us first define the force by dis- 
crete differentiation according to 

F(R)= V&)-V<to , A ( 10 ) 

This definition carries an error O ((A/_R) 2 ). 

In order to reduce this discretization error, 
R should be properly placed within the interval 
R2Y We have exploited the parametrization, 
Eqs. (@,[|!|), to compute this tangential point and 
estimate its error from the uncertainties of the 
fit parameters I = 7/e and f/e. We assumed 
systematic errors of 20% and 100% on these pa- 
rameters, respectively. The discretization effects, 
estimated in this manner, clearly dominate the 
errors at small R. For the subsequent analysis, 
we have chosen 1 <| R2 — R\ |< 2. 

The parametrization dependence of F is dis- 
played in Fig. |^. The circles therein refer to 
Eq. @ and squares to the Cornell form (Eq. (|l3)) 
of the potential. They differ by much less than 
the systematic errors we state. Note that, for the 
purpose of this figure, the (horizontal) shift in the 
tangent point R has been converted into a (verti- 
cal) shift in F. 

3. TWO-LOOP ANALYSIS 

Given F, we base our a g g-analysis on the two- 
loop prediction 



Figure 3. The plateau in A R (R) at (3 = 6.4. The 
inner (outer) part of the errorbars refers to the 
statistical (systematic) uncertainties. 



with t = Iti(ArR). In this part of the analysis, we 
exclude those pieces of data on the force F that 
involve V(R = 1), since the unfolding of lattice 
artifacts is not successful on this point, according 
to Fig. I 

We use Eq. ( [ll|) to convert the force data, 
{F(R)}, into a set {Ap>(R)}, in order to arrive at 
a sensitive test of the two-loop approximation. A 
universal plateau in A#(i?) is found in the range 
Ra < 1 GeV" 1 for all f3 values beyond 6.0. The 
situation for (3 = 6.4 is visualized in Fig. where 
the plateau is found to occur for R values smaller 
than 3.7. The errors refer both to the statisti- 
cal (inner error bars) and systematic uncertain- 
ties (outer error bars). From the data at (3 = 6.4 
and 6.8, we obtain by averaging 



Ar = .660(40)x/ct . 



(12) 



How well does the force scale? In Fig. ^, we have 
compiled all data sets with (3 > 6.0 into one single 
scaling plot of the reconstructed continuum force 
vs. a. To do that we have used appropriate units 



of the string tension v K, as determined from the 
potential analysis in the "long distance domain" , 
r = Ra > .3 fm. 

It is gratifying to observe that scaling is obeyed 
very nicely. Moreover the perturbative two-loop 
formula describes the i?-dependence of the in- 
terquark force impressively well down to energies 
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Figure 4. Scaling plot of the reconstructed contin- 
uum force. The perturbative region is zoomed in 
the inset. Full curve with errorband corresponds 
to the two-loop formula, with A# from Eq. (12). 



as low as .5 GeV, i.e. very close to the Landau 
pole! We emphasize that the f3 = 6.5 data points 
of UKQCD Jl| are in very good agreement with 
these findings. 

To illustrate that our results are not biased 
by the analysis procedure, we finally compare in 
Fig. | the raw lattice potential data with the inte- 
grated perturbative force from Eq. (jll]). We find 
beautiful agreement and no indication of a bias! 



4. SETTING THE SCALE 

It can be argued, that the string tension is 
too elusive a construct and therefore not the best 
ground for laying the foundation to a scale in lat- 
tice gauge theory. To be precise, charmonium or 
bottomium phenomenology is not really sensitive 
to a linear asymptotic rise of the static quark- 
antiquark potential, which will die away in full 
QCD anyhow. SOMMER proposed to base 
the scale fixation on a characteristic length, ro, 
instead, which is connected to the intermediate 
range of the potential and therefore of relevance 
to the physical spectrum of gg-bound states. This 
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Figure 5. Scaling plot of the lattice potential. 
The solid curve is the integrated perturbative 
force from Fig. 4; the dotted curve is —av{R)/R 
with Ay = 1.53Aij. 



length obeys the dimensionless condition 
RlF(R ) = const. = -1.65 . 



(13) 



The rhs is chosen with an eye on phenomeno- 
logical potentials, to tune the value of the "force 
radius" ro = i?o a close to .5 fm. From a Cornell 
paramctrization to the lattice potential 



V(R) = Vo-e/R + KR 



(14) 



we compute ro and map it onto the string tension 
(using the values of the fit parameters at (3 = 6.2), 
with the result 



1.164r " 1 



(15) 
i 



From phenomenology |17| , we retrieve Tq 
400(15) MeV, which induces a value of the rede- 
fined string tension *ftj = 465(17) MeV or, with 
Eq. ([[I), A^T° = 293(18)(11) MeV. One can 
check the string tension data against the precise 
data from GF11 |Q for the p mass. Apparently, 
we cannot get away with quenched calculations 
alone, as we find the ratio m p /^Ja to depend on 
a, at least in the region f3 < 6.2. Extrapolations 
allow for a continuum ratio in the range 1.7 to 2. 
For this reason, we add an asymmetric systematic 
error on the scale which reflects both, quenching 
effects and scaling violations: 



465±i ' MeV 



(16) 
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or A— — 

MS 

5. DISCUSSION 



293(18)163 MeV 



It is instructive to exhibit the amount of scal- 
ing violations remnant in different improvement 
schemes. This is most apparent when plotting 
A^i from the two-loop approximation vs. a. 
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Figure 6. Comparison of various 'improved' cou- 
pling schemes in terms of their apparent two-loop 
A -1 values, as function of a, in units of the string 
tension. 

In Fig. |[ we compare our present results with 
the estimate from the mean- field improved FNAL 
scheme [|| 



fFNAL^/a) =9o 2 (Un) 



.025 



(17) 



We also included the results of two other improve- 
ment schemes into Fig. |^, denoted by gE, gE 2 - 
These "effective" schemes are construed to force 
the perturbative series and Monte Carlo results of 
the plaquette to coincide with each other in lead- 
ing and next- to-leading order, respectively p^,p[. 
The figure provides evidence that a q q from the 
short range force is closest to the asymptotic 
realm of two-loop perturbation theory! In this 
sense, it looks indeed like a "perfect coupling" ! 

In order to extract an estimate for A^gg in the 
continuum limit from the improved schemes, one 
still has to perform extrapolations (see Fig. 
These are rather delicate by their functional form, 



with unknown a and 1/ In a dependencies in- 
volved g). The QCDTARO group j20) suggested 
to flatten A -1 (a) with brute-force, by changing 
to an "optimized" (yet ad hoc\\) renormalization 
scheme through the replacement 



-2 
#FNAL 



— z — Z I 

3fit — ffFNAL + X 



(18) 



where both, xq and the three-loop coefficient in 
the /3-function, are treated as free fit parame- 
ters. It turns out, however, that xq is unphys- 
ically large. 



Figure 7. Comparison of three lattice predictions 
to ajfg for nf = 0. The upper curve is from the 
present analysis, Eqs. (11,12). A : FNAL-MS- 
scheme; □: result of Ref. [7] with its two-loop 
evolution (dotted line). 



In Fig. 0, we display ctj-^ vs. \nq, with the 
momentum in units of the string tension to avoid 
overall scale errors, in the zero flavour world. The 
upper error band refers to the two-loop evolution 
of our a q q values. Two sets of data are shown in 
comparison: the triangles refer to the estimates 
from Eq. (^) , with a determined from the string 
tension, while the point with error bar is the re- 
sult of method (c) , presented by WOLFF to this 
conference [@. Note from Eq. (§) that_ our cou- 
pling scheme is much closer to the MS scheme 
than theirs, as expressed by the relation M 



a- 



MS 



oicdm + l-256a^. DM 



0(a 3 CDM ). 



(19) 
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6. UNQUENCHING 

Having established the method in the pure 
gauge sector, we can now ask the sneaking ques- 
tion 'what will happen in the real world with dy- 
namical fermions present?' In the absence of full 
fledged QCD simulations, we have to resort to 
sneaky answers, i.e. rough estimates || based on 
early day full QCD results. 

For that purpose, we will exploit the work of 
the MT C collaboration who performed a match- 
ing of the quenched and unquenched (4 flavour 
staggered fermions) potentials |2l| by tuning (3. 
Their results enable us to compute the FNAL-A/5 
couplings, in the zero- and four- flavour worlds. 
From the latter, we can extract the matching 
of the accompanying A- values, and find A^L = 

(.44 ± -10)A^L We assigned a 12% (10%) error 
due to matching ( scaling violation) uncertainties. 
The degeneracy of the lattice quarks leads us to 
increase the systematic error: 

A^=(.44ll°)A^=129±8l«MeV . (20) 

This can be converted into an a-value in the five 
flavour world at the Zq mass: 

<*mbn z ) = .1020t& . (21) 

This lattice estimate appears to be rather low, 
compared to the average value from LEP experi- 
ments, which reads d^- p = .118(7). 

' MS v ' 

7. CONCLUSION 

We have seen that a direct determination of 
the strong coupling from the interquark force is 
feasable with the computing power of an 8K Con- 
nection Machine CM-2. The final number for 
a s (mz), shown in Eq. (pl|), carries errors which 
are comparable to the experimental ones. At this 
stage, by far most of this estimated error is due 
to scale uncertainties (15% in A— —°, mainly from 
quenching) and the cavalier — s- 4 flavour conver- 
sion (30%). 

The method is applicable to full QCD. Dis- 
cretization effects can still be reduced in the 
present framework by going to finer lattices. Fu- 
ture effort should also go into the direction of a 
full lattice perturbative analysis. 
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